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The main objective of this work is to investigate theoretically how tilting of an easy axis of a 
single-molecule magnet (SMM) from the orientation collinear with magnetic moments of the leads 
affects the switching process induced by current flowing through the system. To do this we consider 
a model system that consists of a SMM embedded in the nonmagnetic barrier of a magnetic tunnel 
junction. The anisotropy axis of the SMM forms an arbitrary angle with magnetic moments of the 
leads (the latter ones are assumed to be collinear). The reversal of the SMM's spin takes place due to 
exchange interaction between the molecule and electrons tunneling through the barrier. The current 
flowing through the system as well as the average z-component of the SMM's spin are calculated in 
the second-order perturbation description (Fermi golden rule). 
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C-| , I. INTRODUCTION 

Jo . 

Single- molecule magnets (SMMs) draw attention as potential candidates for devices which can combine conven- 
tional electronics with spintronics^. Characterized by a relatively large energy barrier for the molecule's spin reversal, 
a SMM can be used at low temperatures as a molecular memory cell^. For these reasons, transport through SMMs is of 
C$ ' current interes^"^. It has been shown that the molecule's spin can be reversed by a spin current (also in the absence 
of external magnetic neld)&£. The phenomenon of current-induced spin switching is of great importance for future 
applications. Furthermore, it is now possible to investigate experimentally transport through a single molecule^&i^. 
However, using present-day experimental techniques, one can hardly control the orientation of the molecule's easy 
£h ! axisii. 

The main objective of this paper is to investigate theoretically how tilting of the easy axis of a SMM from the 
i— i orientation collinear with magnetic moments of the electrodes (leads) affects the switching process and current flowing 
■ through the system. 
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The system under consideration consists of a SMM embedded in a nonmagnetic barrier between two ferromagnetic 
, electrodes. Electrons tunneling through the barrier can interact via exchange coupling with the SMM, which may 
result in magnetic switching of the molecule. Furthermore, we assume that the spin number of the SMM does not 
C~-- , change when current flows through the system, i.e. the charge state of the molecule is fixed. In the case considered 
here, the anisotropy axis of the molecule (used as the global quantization axis z) can form an arbitrary angle <f> with 
magnetic moments of the leads. To simplify the following description, we neglect the influence of exchange interaction 
with the leads on the ground state of the molecule. Such an influence, however, can be included via an effective 
KnJ ' exchange field. 

£2 . The full Hamiltonian of the system reads 

_ _ i 

H = Hsmm + + Hr + Ht- (1) 



II. MODEL 
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FIG. 1: (color online) Schematic picture of the system under consideration for two collinear configurations of the leads' magnetic 
moments: parallel (blue arrows) and antiparallel (red arrows). The axis y (y 1 ) is normal to the xz (x'z') plane. 
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The first term describes the free SMM and takes the form 



H 



SMM 



= -DSi 



(2) 



where S z is the z component of the spin operator, and D is the uniaxial anisotropy constant. The next two terms of 
the Hamiltonian Ti correspond to the two ferromagnetic electrodes, 
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for q = L (left) and q — R (right). The electrodes are represented by a band of non- interacting electrons with the 
energy dispersion e^ a , where k denotes a wave vector, a is the electron spin index (a = + for spin majority and 
a = — for spin minority electrons), and a^ a (a^c) ^ s the relevant creation (anihilation) operator. Finally, the last 
term of the Hamiltonian H stands for the tunneling processes and is given by the Appelbaum Hamiltonian 1 - rotated 
by the angle <f> around the axis y' — y (see Fig. [T]), 
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where cr = (a x , a v ,a z ) are the Pauli matrices, S is the SMM's spin operator, and a = —a. The first term in 
the above equation describes exchange interaction of the SMM and electrons in the leads, with J q> qi denoting the 
relevant exchange parameter. For the sake of simplicity, we consider only symmetrical situation, where Jl,l — Jr.b. — 
Jl,r = Jr.l = J- The second term in Eq. ([4J represents direct tunneling between the leads, with Td denoting the 
corresponding tunneling parameter. We also assume that J and Td are independent of energy and polarization of the 
leads. Finally, N q (q = L, R) denote the number of elementary cells in the g-th electrode. 

III. THEORETICAL DESCRIPTION 

The electric current / flowing through the system is determined with the use of the Fermi golden rule^£. Up to the 
leading terms with respect to the coupling constants J and Td, the current can be expressed by the formula 



/= T e2 [l^l' + l^ 032 ^ 2 )] (DfDf + DfDf) V 
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+ ^^(DfDf + DfDf)[C(D(r]2m + l) + rieV)+C(D(ri2m + l)-rieV)] \, (5) 
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where e is the electron charge (for simplicity assumed e > 0, so the current is positive for electrons tunnelling from 
the left to right). In the above equation D q a is the density of states (DOS) at the Fermi level in the g-th electrode 
for spin a, and (S%) — J2 m m2 Pm, where P m denotes the probability of finding the SMM in the spin state \m). The 
voltage V is defined as the difference of the leads' electrochemical potentials, eV = — I^r- Finally, we introduced 
the notation: A±(rri) = 5(5 + 1) - m(m± 1), and £(e) = e[l - exp(-e/3)] _1 with = k B T. 

To compute numerically the current / from Eq. ([5|), one needs to know the probabilities P m . To determine them, 
the SMM's spin is assumed to be saturated in the initial state | — 5), and then voltage growing linearly in time is 
applied 6 . Since the reversal process occurs through all the consecutive intermediate spin states, the probabilities can 
be found by solving the set of relevant master equations, 
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FIG. 2: (color online) The average value of the SMM's spin, {S z ), and the current I flowing through the system as a function 
of the voltage V, calculated for the parallel configuration of the electrodes' magnetic moments with P L — 1 and P R = 0.3. 



for —S<m<S, where c = V/t is the speed at which the voltage is increased. The parameters 7™ describe the 
rates at which the spin z component (m) is increased (decreased) by one. These tunneling rates have been calculated 
from the Fermi golden rule and have the form, 
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IV. NUMERICAL RESULTS AND DISCUSSION 



Numerical results have been obtained for the molecule Fes 1 ^ corresponding to the total spin S — 10, whose 
anisotropy constant is D = 0.292 K. Apart from this, we assume J w T<j w 100 meV. Calculations have been 
performed for the temperature T — 0.01 K, which is below the molecule's blocking temperature Tg = 0.36 K, and 
for c = 10 kV/s. It has been also assumed that both the leads are made of the same metallic material characterized 
by the total DOS D = D q + + D q _ w 0.5 per electron-volt and per elementary cell, where denotes the DOS of 

majority (minority) electrons in the q-ih electrode. Furthermore, the q-th. electrode is described by the polarization 
parameter P q defined as P q = (D+ — D q _)/(D q L + D q _). The following discussion is limited to the case, where one 
electrode (the left one) is fully polarized, P L = 1, whereas the polarization factor of the second electrode can vary 
from P R = (nonmagnetic) to P R = 1 (half-metallic ferromagnet). 

In Fig. 2 we show the average value of the z-component of the SMM's spin, (S z ), and the charge current /, 
calculated for several values of the angle <j> an d for parallel magnetic configuration of the leads. The case of <fi = 
((f) = it) corresponds to the situation when the initial SMM's spin is antiparallel (parallel) to the leads' spin moments. 
One can note that the influence of current on the molecule's spin gradually disappears as the angle <f> approaches 
7r. For cf> < 7r/2, the molecule's spin becomes switched from the state | — S) to the state \S). The switching time, 
however, becomes longer and longer as the angle <j> approaches (f> = n/2. At (f> = n/2, which corresponds to the 
situation with the SMM's easy axis perpendicular to the leads' magnetic moments, different molecular spin states 
|m) become equally probable with increasing voltage, and therefore (S z ) — > 0. This is a consequence of the fact that 
when the voltage exceeds the activation energy for the spin-flip process^, the SMM undergoes transitions to upper 
and lower spin states with equal rates 7+ = 7" (see Eq. 

When <f> > 7r/2, the spin state of the molecule is only weakly modified by current, and remains strictly unchanged 
for <j) — 7r. The absence of switching by positive current at large values of <f> (for the assumed parameters) is consistent 
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FIG. 3: (color online) The average value of the SMM's spin {S z ) (a), and the current / flowing through the system (b) as a 
function of the angle 4> m the parallel (P) (a,b) and antiparallel (AP) (c,d) magnetic configurations for V — 2 mV and P L = 1. 



with the conclusion of Ref . [6] , where for collinear configurations and positive current only switching from | — S) to 
\S) states was allowed, whereas positive current had no influence on the state \S). 

The SMM's spin can be reversed due to exchange interaction with tunneling electrons. The latter flip their spins 
and hence add to or subtract some amount of angular momentum from the molecule. As the angle (f> grows, the spin 
orientation of tunneling electrons 'seen' by the molecule and consequently also the transition rates given by Eq. ([7]) 
change as well. Figure [5Jd shows the current flowing in the system as a function of the bias voltage. This current 
strongly depends on the orientation of the SMM's easy axis. This dependence is a consequence of the fact that 
the dominant contribution to current is due to the exchange term (first term in Eq. ^ ), which is sensitive to the 
orientation of the SMM's spin. The curves for = and (f> = tt overlap (except for a small voltage range where the 
switching for (f> = takes place - not resolved in Fig. 

Figure [3] presents the spin z-component (S z ) and the current / in both configurations of the leads' magnetic 
moments, plotted as a function of the angle cf) and calculated for V — 2 mV. For tfi < ir/2, the spin switching takes 
place in both parallel and antiparallel magnetic configurations. Figure [3] also indicates that the current at <f> — n/2 is 
independent of the magnetic configuration as well as on the polarization parameters of the leads. 

In conclusion, we have shown that tilting the easy axis of a SMM from the collinear orientation relative to the 
leads' magnetic moments has a significant influence on the reversal process of the molecule's spin, as well as on 
current flowing through the system. 

Acknowledgements This work was supported by funds from the Ministry of Science and Higher Education as a 
research project in years 2006-2009. One of us (MM) also acknowledges support from the MAGELMAT network. 



* Electronic address: barnas@amu.edu.pl 

1 D. Gatteschi and R. Sessoli, Angew. Chem. Int. Ed. 42, 268 (2003). 

2 C. Joachim, J. K. Gimzewski and A. Aviram, Nature 408, 541 (2000). 

3 C. Timm and F. Elste, Phys. Rev. B 73, 235304 (2006). 

4 G.-H. Kim and T.-S. Kim, Phys. Rev. Lett. 92, 137203 (2004). 
6 F. Elste and C. Timm , Phys. Rev. B 73, 235305 (2006). 

6 M. Misiorny and J. Ba rnas , Phys. Rev. B 75, 134425 (2007). 

7 F. Elste and C. Timm, |cond-mat/061li08| 

8 H. B. Heersche et al., Phys. Rev. Lett. 96, 206801 (2006). 

9 M.-H. Jo et al., Nano Lett. 6, 201 4 (2006). 

J. J. Hen derson et al., |cond-mat/07030"i"3] 

1 C. Timm, |cond-mat/0702220| 

2 J. Appelbaum, Phys. Rev. 17, 91 (1966). 

3 W. Wernsdorfer and R. Sessoli, Science 284, 133 (1999). 



